Statistics of current activity fluctuations in asymmetric flow with exclusion 
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We consider steady-state current activity statistics for the one-dimensional totally asymmetric 
simple exclusion process (TASEP). With the help of the known operator algebra (for general open 
boundary conditions), as well as general probabilistic concepts (for the periodic case), we derive 
and evaluate closed-form expressions for the lowest three moments of the probability distribution 
function. These are confirmed, to excellent degree of accuracy, by numerical simulations. Further 
exact expressions and asymptotic approximations are provided for probability distributions and 
generating functions. 
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I. INTRODUCTION 

The one-dimensional totally asymmetric simple exclu- 
sion process (TASEP) is a biased diffusion process for 
particles with hard-core repulsion This model 

exhibits many non-trivial properties, and is considered 
a paradigm in the field of non-equilibrium phenomena. 
While a large number of exact results are known for the 
steady-state, ensemble-averaged, values of relevant quan- 
tities, such as the systemwide current and global density, 
the study of the corresponding fluctuations has proved 
to be rather more intricate. To quote only a few relevant 
developments, exact expressions for the diffusion con- 
stant were found for systems with periodic (PBC) (y| and 
open [t'I boundary conditions (BC); the full probability 
distribution function (PDF) of current fluctuations was 
similarly considered for both PBC [8] and open |9] BC. 
Very recently, a number of new results have been found 
for current fluctuations in systems with open BC p^l - [l3| . 

In this paper we investigate steady-state current activ- 
ity fluctuations in the one-dimensional TASEP for both 
periodic and open boundary conditions. It is important 
to recall at the outset that the current activity is not iden- 
tical to the standard current (although the first moments 
of the respective distributions coincide) . As explained in 
detail below, the former quantity is static, in this sense 
akin to the instantaneous (local or global) particle den- 
sity, while the latter is a dynamic one. Our results ex- 
tend and complement earlier analytic work on the joint 
current-density distribution for the TASEP [T3 . [l5 | . 

For the problem of flow with exclusion, the time evolu- 
tion of the 1 -t- 1 dimensional TASEP is the fundamental 
discrete model. The particle number ni at lattice site ^ 
can be or 1, and the forward hopping of particles is 
only to an empty adjacent site. The current across the 



bond from £ to £ + 1 depends also on the stochastic at- 
tempt rate, pi, associated with it and is thus given by 
Jt,t+i = Pene{l — n^+i) . For the homogeneous case of 
P{ = p which is considered here, one can effectively make 
p = 1 in numerical simulations, provided that the inher- 
ent stochasticity of the process is kept, via, for example, 
random selection of site occupation update [l^. This 
amounts to a trivial renormalization of the time scale, 
and is the procedure followed in all numerical work re- 
ported in this paper. 

For PBC the main parameter is the fixed density, which 
is uniform in the steady state. With open boundary con- 
ditions, one has the following externally-imposed param- 
eters: the injection (attempt) rate a at the left end, and 
the ejection rate /? at the right one. The phase diagram 
in a-/3 parameter space is known exactly, as well as many 
other steady state properties [iHi, [13, [ll, [O, [H] . 

In the open-boundary case both current and global 
(position-averaged) particle density are fluctuating quan- 
tities (as opposed to the case with PBC for which the 
density is fixed). On the other hand, while continuity dic- 
tates that the average stationary current on every bond 
must be the same, steady-state density profiles may be 
non- uniform fs]. 

For PBC the total (instantaneous) activity A within 
the system is defined as the number of bonds that can 
facilitate a transition of a particle in the immediate fu- 
ture. Thus it equals the number of pairs of neighboring 
sites that have a particle to the left and a hole to the 
right [3, Foi' systems with open BC, an alternative 
definition includes also the injection and ejection bonds 
at the system's ends, though these have to be weighted 
by the respective injection and ejection rates, a and /3. In 
the following, for open BC we always include the contri- 
butions given by the latter bonds, unless explicitly stated 
otherwise. 

Thus for an L-site system with PBC, one has: 
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^ = ^7i,(l-"£+i) (PBC) 
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(with riL+i = ni), while for the open BC case (L sites 
and L + I bonds, including the injection and ejection 
ones), one has: 

L-l 

A = a (1 - + ^ rt£ (1 - Ui+i) + PriL (Open BC) . 
1=1 

(2) 

The activity is, therefore, a snapshot of the system at a 
given moment in its evolution; in this sense, it is as much 
of a static quantity as, for example, the instantaneous 
global density. By contrast, the current is a dynamic 
object, as it reflects the stochastically-determined parti- 
cle displacements which actually take place during a unit 
time interval. 

The detailed balance measure used here for the PBC 
case, and the operator algebra used for open BC, apply 
to the steady state [6[ [3, thus they naturally 
yield the activity; although they can be extended to deal 
with static measures of dynamic quantities such as the 
diffusion constant P, 't'I and more general aspects 0, [l^ > 
we do not attempt similar developments here. 

The qualitative distinction between activity and cur- 
rent has the consequence, upon the (properly normalized) 
moments of the PDFs of the two quantities, that while 
their first moments (averages) coincide, higher cumulants 
differ. For A as defined in Eqs. ([ij, ([2]) one has, respec- 
tively: 

J=Ua) (PBC) ; (3) 



J = {A) (Open BC) , (4) 

where J is the average steady-state current through any 
bond in the system. Equality of averages can be under- 
stood by recalling that successive steady-state snapshots 
(activity configurations) are generated via the interven- 
ing particle hoppings, which constitute realizations of the 
system's current. In our simulations we have verified that 
this property holds, to within numerical accuracy, in all 
cases investigated here. However, the connection at this 
level is not sufficient to warrant equality of higher mo- 
ments of the PDFs. 

Here we consider the fluctuations of the global, in- 
stantaneous activity IT5l |. as opposed to the local, 
integrated current which is the subject of many extant 
works 0-[ll- 

Section [H] below gives the theoretical approach used 
in this work. It first addresses open boundary cases with 
a-|-/3 = 1 and then general a, /3, and then treats PBC. In 
Section Hm the results of numerical simulations are given, 
first for open BC and then for PBC. In Sec. lIVI we provide 
a global discussion of the second and third cumulants of 
the activity PDFs, everywhere on the a — /3 plane (for 
open BC); finally, concluding remarks are made. 



II. THEORY 

A. Open boundary conditions with a + /3 — 1 

For open boundary conditions, steady state configura- 
tions are typically highly correlated. So when (as in the 
pioneering work of Ref. y) such configurations are writ- 
ten as strings of variables D and E {D, E representing 
particle, vacancy respectively) these variables have to be 
operators, with an algebra {D+E = DE) consistent with 
the steady-state properties. In the special case a + (3 = 1 
the correlations disappear Q (configuration probabilities 
factorize, corresponding to a product measure) and D, E 
can be taken as c-numbers d, e respectively equal to the 
steady state averages (n^), (1 — n^). Consequently, along 
the a + f3 = l line these averages are independent of ^ Q ; 
consistency with probability normalization {d + e = 1) 
and current conservation {ea = de = dj3) thus implies 
e = d = a. 

The factorization of probabilities for variables on 
different sites makes this case the easiest to begin 
with. However, different bonds share a site variable if 
they are adjacent, and that makes the calculation of ac- 
tivity statistics beyond the first moment non trivial even 
in this case. 

The average activity {A) is unaffected by the shared 
variable on adjacent bonds, since it is the average of a 
sum of terms each of which contains independent site 
variables and so has a factorizing average. So in this 
case, from Eq. ^ one has: 

{A) = {L + l)aP (a + ^ = 1) . (5) 

Already in the second moment, (A^), the effect of cor- 
relations are present, making it different from {A)'^ . To 
obtain any moment, and also the probability distribution, 
we consider the generating function 

Logarithmic derivatives of (e'^'^) give the cumulants 
Cn of A. The probabilities P{Ai) of the discrete possi- 
ble outcomes Ai for A are the coefficients in the terms 
proportional to e"^"** in the expansion of the generating 
function. 

The average in Eq. ^ does not factorize because of 
the correlation of adjacent bond variables. The follow- 
ing method, involving a transfer matrix along the chain, 
easily handles that. 

Defining Q = 1 — ni, since each Q takes values (0, 1) 
with respective probabilities (a, /3), the average in Eq. ^ 
can be obtained by summing over the possible configu- 
rations of each variable in turn, starting, say, with I — 1 
and working along the chain to £ = L. Using sums 
over the Q variables, with their weights p{Ci) {= ct or 
/3), the first step involves E(^i=o,ip(Ci) e'^"'^' (^"^^^ = 
Q, gA C2 _|_ ^ gAa ^ 'pjiig result can be written in the linear 
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form (62 + C2 C2) where 62 = a + fie^" and C2 = a (e"^ — 1). 
The sums over the next variables then produce similar 
forms, where the relationship between successive ones is 
provided by 

bi+i + a+iCi+i^ J2 PiCi)^^'-'''^'^'^'*'ib, + c,Ci) ■ (7) 

This is the same as the action of the transfer matrix 

'a + (3 P 



a7 ; ' 



(8) 



on the vector ( "'^ ) , where ^ — — 1. 



The consequence is 



(9) 



where, from e"''^ = bi + ci (i, bi = 1, ci = e^" — 1. 

This is readily evaluated using the representation in 
which T is diagonal, with the result: 



e^^) ^a+fj,^ ^ +a-fj, 



.L-l 



where 



1± ^/(a- /3)2 +4a^e-^ 



(10) 



[1±^] (11) 



are the two eigenvalues of T, and the coefficients a+, 
a_ are also known functions of a, /3, 7 (given in Ap- 
pendix 121 • 

This result provides the nth cumulant C„ of A via 



ln(e 



\A\ 



A=0 



From Eqs. (fTU)l and (fTTjl . for large L one has: 



ln( 



L ln/i+ + 0(1) 



(12) 



(13) 



Consequently the asymptotic large- L values of the first 
few cumulants are: 



Ci = (A)«L^ln/i+ 



= La(3 + 0{l)] (14) 



The special case a = /3 = 1/2 is simple and instructive. 
There the above leading large-L results for C3 and skew- 
ness reduce to zero, which turns out to be exactly true 
for any L: for a = /? = 1/2 the coefficient a_ vanishes 
and the full result from Eqs. ([TU| and ([TT|) (exact for any 
L) is: 



1 



;(1 



L-1 



(17) 



From this it is easy to verify that C3 = S* = 0, and that 
Ci = (L + l)/4, C2 = (L - 1)/16. 

The corresponding general L exact results obtained 
from Eqs. (fTU| and (fTT|) for any a, j3 with a -|- /3 = 1 
are given in Appendix together with validity condi- 
tions for the large-L results in Eqs. P^ - p^ . 

For a = /? = 1/2 the generating function result, 
Eq. (|17p can be expanded in the form 



= 2-(^-i)^i_iC;.e^ 



(n+l) 



(18) 



n=0 



corresponding to a binomial distribution of probabilities 
P{Am = m/2) = 2^(^^^) L-iCm-i for the possible out- 
comes m/2, l<m<L, of the activity. 

For general a, /? (still with a + (5 — 1) the possible 
outcomes for A will obviously be linear combinations of 
a, /3, and integers. Their distribution, resulting from 
Eqs. (|10p and (|TT|l . has (for a 7^ /3, where a_ is non- 
vanishing) a "discreteness alternation". This is a direct 
consequence of Eqs. ([TU)) and PT|) . as can be seen from 
the following alternative form: 



L-l 



(19) 



n=0 



A=0 



We conclude this part on activity statistics for a -t- 
/3 = 1 by briefly commenting on local fluctuations. 
The product measure makes these trivial for this case. 
So the injection current activity (from aQi) is a or 0, 
with probabilities /?, a respectively. The local activ- 
ity for the ejection current [/3(1 — Cl)] and for the 
internal current [ (1 — Q) Q+i — ] across any 

specifled bond similarly have binary distributions, with 
(outcomes) (respective probabilities) being (/3,0) (a,/?) 
and (1,0) (a/3,1 — a/3), respectively. So for example 
{Aij+i) = a/3, = = a/3^, indepen- 

dent of and the consequent skewness is (/3— a)(a/3)~^/2. 



C2 = - w L a/3 [1 - 3a/3] + 0{\) ; (15) 

C3 = (A3)-3(^)(^2^-|-2(^)3« 

«La/3 [l-9a/3-|-20(a/3)2] +0(1) .(16) 

Thus the skewness ^ S = C^i/Cl'^ is of order L^^/^ 
for large L. 



B. Open boundary conditions with general a, /? 

Moments of the eurrent aetivity 

For general a, /3 one must make full use of the steady 
state operator algebra to take care of the absence of 
product measure. In this subsection we derive expres- 
sions for low moments of the activity. 
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We exploit the representation of configurations, and 
their probabiUties, using strings of operators D and E, 
representing respectively particle or vacancy at a site. 
Configurations are specified using the variables Q intro- 
duced in Section III Al The probability of the configu- 
ration {Q} for L sites is given, through the associated 
operator string S{{Q}) = l\j=i ((1 " Q) D + Q E), as 

{w\si{c,})\v) 



PiKi}) 



with Zl = (VK|C^ \V) 



(20) 

Here, C = D + E represents both possibilities (occupied 
or vacant) at a given site; and the vectors {W\ and \V) 
are such that 



{W\aE = {W\ , [iD\V) = \V) 



(21) 



Eqs. (pTjl follow from the balance between the current 
A = DE across any internal bond and the injection or 
ejection current (at the left or right boundary), using the 
operator algebra steady-state relation K — C . Using 
Eq. (12) , the nth moment of the total activity is given by 



{0} 



L-1 

E 

1=1 



{1-0.)Q+i+P{1-Ql)T X 



^ {WmUil^Q)D + QE)\V) ^ ^^^^ 

Consider, for example, the simple case n = \. There 
the contribution from (1 — C,i)C,i^i is non-zero only for 
1 — Q = 1, Q+i ~ 1 whose associated probability is 
(using A = DE^C) 



(23) 



Zl Zl 

Correspondingly, using Eq. (|2ip the non-zero contribu- 
tions a from a(i and /3 from /?(1 — (^l) occur with re- 
spective probabilities a^^ Zl^i/Zl and /3^^ Zl^i/Zl, 
making 



{A)^{L + 1) 



Zl-1 



(24) 



This is the standard result generalizing Eq. for a + 
13 = 1, where Zl = {a /Sy^- 

As in Section III Al complications arise in higher mo- 
ments from the shared variable in adjacent bonds. For 
a + /3 = I the transfer-matrix approach overcame these. 
That method is ruled out by the non-commuting vari- 
ables in the present general case. Instead the following 
direct method can be used for low n. 

For n = 2, Eq. ((22|) includes the product of bond 



terms [T,i=ii^-Ce)Ci+i) - Cm)Cm+i) • For 

m — £ 7^ 0,1 in the double sum the sites involved are 
distinct, which simplifies the averages occurring. For ex- 
ample, for m > £ + 1 the average is 



{W\C^-^DE C 



m~£-2 



Zr- 



L-2 



Zl 



Zl 



(25) 



This result applies for any £, m satisfying m > + 1, or 
£ > m + l (altogether (L — 2)(i — 3) terms of the double 
sum). For each of the L—1 terms with £ = m the product 
of bond variables is ((1 — C^)C^_|_i) = (1 — 0)0+i, having 
average {W\C^'~^ A C^-(^+i) = Zl-i/Zl- 

For the 2(L — 2) terms with £ = m ± 1 the product 
of variables includes both a 1 — C and a C for the shared 
site, whose product is zero. Hence, for the part A' = 

J2e=ii^ ^ 0)C^+i of the activity coming from just the 
internal bonds, 



{A'') = (L'-5L + 6)^ + {L-l)^ . (26) 



Including also the injection and ejection contributions 
gives 



{A^) = {{aa+A' + (3il-CLW = 

{a^Ci + 2a/?Ci(l - Cl) + 13^1 - Cl)> + 
+ {2 (aCi + /3il - Cl)) A') = 
Zl-1 



Zl 



2T 



(27) 



where 



L-l 



r = ({aCi(i-Ci)C2 + «Ci E(i-C^)C^+i 



e=2 



L-2 



(1 - Cl-i)Cl /3(1 - Cl) + /3(1 - Cl) J2^l - 0)0+i}> 



2(L-2)^. (28) 



The case n = 3 is in principle a straightforward gen- 
eralization of the above. However, the {A'^) part of it is 
considerably more complicated. Details of the evaluation 
are given in Appendix |B] 

The collected results for the first three moments of 
both A and A' are: 



(A) ={L + 1) 
{A') = iL-l) 



Zl- 



1 



Zl 
Zl-1 



(29) 



n = 2 : 

(A2)^L(L-l)^ + (L-l + « + /?)^ 

^L ^L 

(A'^) = (L - 2)(L - 3) ^ + (L - 1) ^ 



(30) 
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n = 3 : 

{A') = {L-l){L~2){L-3)^ 



Zr 



+3{L - 1){L - 2 + a + 13) 



Zl- 



Zl 



{A") = iL-3)iL-A){L-5) 



Zl 

Zl-3 

~zr 



+3(L-2)(L-3)%i + (i-l)%^ . (31) 

The functions Zl are available from Ref. [31, or by using 
the generating function method of Refs. [ij andfla. The 
general expression is: 



1. LiiL-ey. ^ 



(32) 



fe=0 



Eq. p2p reduces to simple closed form expressions in 
the cases: 



(i) a + /3 = 1 : 

(ii) a = ^ = 1 



Zl- 
Zl 



(2L + 2)! 



(33) 
(34) 



(L + 2)!(L + 1)! ■ 

From Eqs. and one can show that, for a = 

/3 = 1 and large i, the skewness S (of the distribution 
for A) varies as 



S^<iL-'^''^{l + 0{L-^)) (a = /3=l) 



(35) 



In Section HITB] below we describe numerical simu- 
lations performed for the cases covered by Eqs. ((33)) 
and p4p . We return to the general a, (3 case, using 
Eq. dnil), in Section [IVl 



2. The activity distribution - I 

This section is concerned with the internal bond ac- 
tivity distribution, for general a, (3. We again use the 
internal bond current activity operator A — DE, and 
C = D + E, which represents both possible single-site 
configurations. 

For n successive sites, C" represents all possible con- 
figurations; we generalize this by representing by (C") 
the subset which contains no A's. Hereafter in this 
subsection the round bracket pair ( ), when enclos- 
ing a string of C"s, or a linear combination of such 
strings, will represent the formal property of picking 
out the subset containing no A's. Then, by defi- 
nition, (C"'« ) A (C"! ) A (C"-2 ) A • • • (C"""-i ) A (C"~" ) con- 
tains many configurations, all having n active internal 
bonds (total internal activity A' = n), from the n A's in 
a run of L = 2n + X]m=o sites. Then, 



SL,A'{{Cn}^Y.--- E ■■■ (C'^--'-)A(C^-') 

(36) 



ro>0 r.,>0 



with To -I- ri + • • • + TA' = L — 2A', contains all configura- 
tions with internal activity A' in an L-site system. So, in 
such a system the probability of internal current activity 
A' is 



PL{A') = zj:'{w\SL.A'{{cn}\v) 



(37) 



From their definitions it follows that (C") and Sn,A' are 
related by 



(C") = C" - £ Sn,A' , 



(38) 



where [X] denotes the integer part of X . The steady 
state algebra relation A = C then makes (C") a function 
of C only, so no lack-of-commutation difficulties arise 
and one finds that Sl,A' is the coefficient of 7 



L-2A' 



the expansion of 



A' + l 



That makes Sl = 



12 A' =1 ^L.A' the coefficient of 7-^ in S'(7), where 



S{l) = 



1 



1 -7C 



1-72C 



1 



1 -7C 



(39) 

Then, consideration of Y1^=q 7™ eventually gives the 
following reduction of the formal operation defined by 

( ): 



(40) 



,l-7Cy l-C(7-72) 
Hence, Sl.A' is the coefficient of 7-^"^^ ir 

5'(7,A')=C^' [l-C(7-7')]^*^" 



(41) 



The required coefficient can be found by the expansion 
of the function of C in powers of (7 — 7^), and then using 
the binomial expansion of each power of (1 — 7) occurring. 
Taking the matrix element {W\ ■ ■ ■ \V) of the resulting 
C-dependent coefficient and using the definition of Z^ 
then provides the following exact result for the internal 
activity distribution: 



Pl{A') 



1 



E 

s=0 



Zl-A'-s{-IY 



{L-A'- sy. 



A'\s\{L-2A' -2sy. 

(42) 

An alternative expansion strategy from Eq. (|41l) onward 
is available for the special case a + (3 — 1. Because the 
operators reduce to c-numbers in this case, C in that 
Equation becomes (a/3)^^, so [1 — C{j — 7^)] becomes 
(a/3)-i(7-7+)(7-7_) where 7± = i(l±Vl-4a/3) = 
a, /3 (the same as the eigenvalues /i±(A) of the transfer 
matrix of Section IH Al at A = 0) . 

In particular. 7+ = 7- = 5 for a — (3 
that 



i , so it follows 



Pl{A') = 2-'^ L + lC2A' + l 



(43) 
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For large L, where boundary current contributions to A 
are less by a factor 0{1/L) from the total internal bond 
current contributions which comprise A', Eq. (|43p be- 
comes the same as the result from Eq. p^ . 

Though the method just presented applies for any a, 
/?, it is not suitable for inclusion of boundary current 
contributions. We briefly proceed next to a more so- 
phisticated approach to the current activity distribution 
which is in no such way limited, and is related to pre- 
vious work [fl on the joint distribution for number 
and current activity. 

3. The activity distribution - II; Generating functions 

The method of Refs. [ij and fisl represented any micro- 
scopic configuration of N particles with internal activity 
A' in the L-site open boundary TASEP as a sequence 
of A' objects of form DP' E''^ pj, hj > 1, with ho > 
-E's to the left and pa > D^s to the right because of 
the possible injection and ejection. The operator string 
representation for the set of all such configurations is ob- 
tained by summing over all {hj, pj} as specified above, 
subject to the additional constraints coming from speci- 
fied L, iV, and A': 

A' A' 

Y^{h,+p,) = L, Y.Pi^^- (44) 

/iQ > (or po > 0) corresponds to unit injection (or 
ejection) current activity, accounting for any difference 
between A' and A. The steady state measure is then 
obtained from the {W\ ■ ■ ■ \V) matrix element of the op- 
erator string. This was reduced in Refs. [13 and [l5l using 
a generalized version of the original operator algebra of 
Ref. ^. 

We are now interested in just the activity distribu- 
tion, for specific i, so the constraint on N drops out 
and the remaining constraints can be enforced using an 
integral representation of a (5-function. The resulting ex- 
pression for the total activity probability function Pl {A) 
can then be reduced using the modified operator alge- 
bra of Refs. [13 and[il The result is that Pl{.A) is the 
coefficient of 2;^~2-4 jj-^ 

^("'^'-^'^)-^ (z-a)(z^"/3)(l-2z)2^ ' 

where, in the numerator, a' = {a — z)/{\ — 2z), /?' = 
(/3 — z)/{l — 2z) replace a, (3 to obtain the "renormal- 
ized" version from the usual Zl = ■ Equivalently, 
the double-integral representation of Refs. [13 and[l5| for 
the joint distribution Pl{A,N) can be summed over N 
to provide the above form. 

As well as giving rise to the "renormalized" Z'^ ^ , the 
modified algebra also gives a generating functional for 
Z2^', that can be exploited following Ref. [H to give a 



direct form for the generating function F" {X) = {e^^) 
as the coefficient of z^~^^ in 

Uf{z, A) = ^ [{2a - 1 + .^(l-2z)2-4z2eA} x 

X {2/3 - 1 + y^(l-2z)2-4z2 e^}]-^ . (46) 

As expected, for a + (3 = I the square bracket simplifies, 
and in the special case a + (3 = ^ it becomes a product 
of two factors linear in z. From this it is easy to recover 
Eq. dni) in Section |lT3 

The form obtained from Eq. (j45p for the distribution 
function is not easy to simplify, except in the asymptotic 
large L limit. There, saddle point integrals analogous 
to those in Refs. [13 and [l^ give the activity distribu- 
tion around its peak, apart from A-independent factors 
or others subdominant at large L. 

The asymptotic results for all phases can alternatively 
be obtained by integrating out the density from the 
asymptotic joint distribution results obtained in Refs. [13 
and 15. So, the activity distribution can be written for 
large L in the form 

Pl(A)= f dpexp[~Lg{p,j)] , (47) 
Jo 

where for the low current phase with a < /3 
g{p, j) = 2j In j + ip- j) \n{p - j) + 
+ {1- p-j)ln{l- p-j) + 
+pln((l-a)/p) + (l-p)ln(a/(l-p)) , (48) 

and for the maximal current phase g{p,j) is of similar 
form but without the last two terms in Eq. (j48p . 

The leading large-L approximation to Eq. P7)l . pro- 
vided by the Laplace method (i. e., exp[— L g(p, j) ] where 
p(j) in [0, 1] minimizes g) gives, for example, a maximal 
current phase PDF independent of a and /3, and the same 
as the Laplace result for the low current phase PDF at 
a = P = ^, which is consistent there with the general-L 
result under Eq. ([TS]). 

Though of theoretical interest, the results for probabil- 
ity distributions obtained by the methods of the last two 
sections are of less value for detailed comparisons with 
numerical simulation results than the exact tractable ex- 
pressions provided, for any L, in Section IIIB II for low 
moments of the activity. See also Section IIVI 

C. Periodic boundary conditions 

For PBC, the general ideas concerning operator 
strings, recalled in section IIIB 1[ are again applicable. 
A number of simplifying features occur, besides having a 
fixed number of particles in the system. From the very 
simple form of detailed balance equations applying in this 
case [6] , for a given number M of particles all steady-state 
configurations C are equiprobable: 

P{C) = (PBC) . (49) 
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We introduce the quantity 



III. NUMERICAL RESULTS 



_ [M -n){L-M -n) 
an = ; : 



(50) 



for n = 0, 1, 2, which wih prove convenient in what 
fohows. 

With the notation estabhshed in Section fll Al the av- 
erage of the two-site function (1 — ) Q+i is independent 
of i, and proportional to the number of combinations in 
which site £ is occupied, site ^ -I- 1 is vacant, and the 
remaining Af — 1 particles are distributed among the re- 
maining L — 2 sites. So, 



{A)=L{{l-Ci)Ci+i) 



lCm 



(51) 



For (A^), similarly to the discussion in Section fllB 11 
one needs averages of four-operator products, namely 
Xtm = ((1 - 0) Ct+i (1 - Cm) Cm+i>- Three possibilities 
arise: 

(i) m £ [L terms] : Xu = ((1^0) Ci+i) = ".q/L [using 
= Ci, etc, and Eg. ([ST])]; 

(ii) m = i±l [2L terms, from the ring geometry ] : 
X^i±i = [since (1 - Q) 0=0]; and 

(iii) m £, £ ± 1 [L^ —3L terms]: contributions come 
from having sites £, m occupied, £-\-l, m + \ vacant, and 
M — 2 particles distributed among the remaining L — 4 

sites, so Xtra =L-4 Cm-2/lCm- 

Adding up (i)-(iii): 

(^2) =ao(l+ai) . (52) 

Details of the evaluation of {A^) are given in AppendixICl 
The result is: 



{A^) =ao(l + 3ai+aia2) . 



(53) 



Introducing the variables p = M / L\ p = p{l — p); A = 
pL'^; and L = L—1, Eqs. ((5T|l - (|53)l give the cumulants of 
the PDF in a form which evinces the expected particle- 
hole duality (invariance under p O 1 — p): 



L 



Co = 



A [A - L] 



LHL-l) 

A [4A2 - A(£2 + 6L) +L^{L + 2)] 
L3 (Z-l)(L-2) 



(54) 



(55) 



(56) 



For large L. the skewness S is given, to leading order in 
L-\ by 



5 = /i-i(4^-l)L-i/2 



(57) 



provided that p ^ 1/4, i.e. p ^ 1/2. In the latter case, 
consideration of higher-order terms shows that 



{p = 1/2) . (58) 



A. Introduction 

We considered lattices with L = 2™ sites. We found 
that the size range corresponding to 5 < to < 8 is gen- 
erally suitable to highlight L-dependent effects which 
vanish for L — >■ 00, and for which we can provide ac- 
curate numerical checks of the theory developed in Sec- 
tions HEl [nBll and HTCl On the other hand, such sizes 
are large enough to prevent lattice discreteness from play- 
ing a significant role [ except for the " discreteness alterna- 
tion" effect referred to above, in connection with Eq. ([T^ . 
which for systems with open BC is a general feature un- 
less a — P — ^ ov \ \. For very low densities with PBC, 
for reasons explained in Section llll Cl below. we extended 
the size range up to m = 11. 

A time step is defined as a set of L sequential update 
attempts, each of these according to the following rules: 
(1) select a site at random; (2a) if the chosen site is the 
rightmost one and is occupied, then (3a) eject the particle 
from it with probability /3; alternatively, (2b) if the site is 
the leftmost one and is empty, then (3b) inject a particle 
onto it with probability a; finally, if neither (2a) nor (2b) 
is true, (2c) if the site is occupied and its neighbor to the 
right is empty, then (3c) move the particle. 

Thus, in the course of one time step, some sites may 
be selected more than once for examination, and some 
may not be examined at all. 

For the various sets of a, (3 (for open BC) or values of 
p (for PBC) considered here [because of particle-hole du- 
ality, we take only a < j3 for the former case, and density 
p < 1 /2 for the latter [ , we ascertained that , starting from 
an initial random configuration of occupied sites (usually 
with overall density p = 1/2), riin = 4, 000 time steps was 
long enough for steady-state flow to be fully established, 
for all system sizes L < 256. This also applies to the 
special case of systems with PBC, L < 2048 and very 
low densities, examined in Section IlII CI 

We collected steady-state current activity data (typi- 
cally iVsam = 10^-10^ independent samples) to produce 
the respective PDFs. Accurate evaluation of their mo- 
ments of n-th order, 1 < ?t. < 3, involves running A'set 
independent sets of A'sam samples each; from the spread 
among the averaged moments for the distinct sets, one 
then estimates the root-mean-square (RMS) deviation of 
each relevant quantity. As is well known [21j[, such RMS 
deviations are essentially independent of iVget as long as 

— 1/2 

A'sct is not too small, and vary as A^sam • We generally 
took Nset = 10. 

Our results for the activity PDF are given in terms 
of the reduced variable x = {A — {A))/Lq [Lq = L + 1 
for open BC, or L for PBC[ which, as a consequence of 
Eqs. ^ and Q, is suitable to highlight similarities and 
differences between the fluctuation statistics of activity 
and those of the standard current J . 
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-0.05 0.05 
X 

Figure 1. (Color online) PDF for the reduced variable x = 
{A-{A))/{L + 1) for a = /3 = 1/2, L = 256 [see Eqs. (gj, ©]. 
The points are from numerical simulations, A'sam = 10*' sam- 
ples. The full line is the binomial distribution given following 
Eq. (UHll. 

B. Open boundary conditions 

Initially we consider the a + /3 = 1 line. For the spe- 
cial point a — fi = 1/2, the considerations of Sec. Ill Al 
[see also Eq. (gS]) in Sec lIIB2l ] indicate that the PDF 
is a binomial function. Fig. [T] shows the corresponding 
numerical simulation data for a system with L = 256 
sites, together with the appropriate binomial distribu- 
tion (which, for this value of L, is indistinguishable from 
a Gaussian function of a continuous variable). The fit 
between the two sets is indeed excellent. 

On the a + f3 — 1 line, the discreteness alternation ef- 
fect associated with Eq. is illustrated for a = 1/4, 
/3 = 3/4 in Fig. [21 To avoid smoothing out of the effect, 
binning of the variable x = {A ~ {A))/{L + 1) must be 
carefully chosen; for the simulations reported in Fig. [5] we 
used bins of width [2(L -f 1)]^^. In order to give higher 
resolution to the central part of the distribution, where 
the sawtooth pattern is quantitatively more significant, a 
linear scale has been used on the vertical axis of the Fig- 
ure. The Gaussian fit shown in Fig. [5] takes into account 
the full set of data, thus averaging over the oscillations. 
Remarkably, it produces a rather accurate fitted value, 
A = 0.0178(3) for the corresponding width of the distri- 
bution (see Table H] below for comparison). 

Still on the a + j3 — \ line, for small a and large /3 
(or, by particle- hole duality, for large a and small /3), 
Eqs. (fT5)l and (fTO)) show that for large L the skewness 
S is positive and large, proportional to (a/?)^^/^. The 
average activity being {A) = [L -\- 1)q;/3, the minimum 




-0.05 0.05 
X 

Figure 2. (Color online) PDF for the reduced variable x = 
{A - {A))/{L + I) for a = 1/4, /3 = 3/4, L = 256 [see 
Eqs. Q, ((5l)[. The points are from numerical simulations, 
A'^sam ~ 10^ samples, and exhibit the discreteness alternation 
implied by Eq. p9|l . The full line is a Gaussian fit to the full 
set of simulation data. 



allowed value of the quantity x ~ {A — {A))/{L + 1) is 
Xmin ~ — a/?! i-e., |a;min| <C 1 for a <C 1, while positive 
values of X < 0.5 are permitted. This is illustrated in 
Fig. [21 where the discreteness alternation has been av- 
eraged out by using bins of width {L + 1)^^ for the x 
variable. The calculated values of S for the curves shown 
are = 0.4471 . . . , 0.3158 . . . , and 0.2232 . . . , respectively 
for L = 64, 128, 256 [see Eqs. dM])-® ]■ Large, pos- 
itive skewness is expected to be a general feature also 
away from the a + jS = 1 line, anywhere for small a, 
large /? (and the complementary region at large a, small 

/?)■ 

Away from the a + /3 = 1 line, we report data for 
a — (3 — \, where the averaged moments are given by 
relatively simple expressions [see Eq. ([M| ]. The activity 
PDF is shown in Figure [H together with the result of a 
Gaussian fit to the data. Visual inspection indicates that 
the fit is of a similar quality to that for the a = /3 = 1/2 
data, displayed in Fig.[T] Indeed, according to Eqs. (|^^ - 
(l?T|) . ([M| . and although a non- vanishing skewness 
must be present in this case, it is quantitatively rather 
small for large L [see also Table 

Results from numerical evaluation of the first three mo- 
ments of the activity PDF, for the values of a, /3 discussed 
above, are given in Table HI together with the theoretical 
predictions summarized in Eqs. The agree- 

ment is very good, except that accuracy for the skew- 
ness is comparatively low. This is especially true for 
a = /3 = 1 , in which case much narrower error bars would 
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Figure 3. (Color online) PDFs for the reduced variable x = 
{A — (A))/ (L+ 1) for systems with open boundaries, a — 0.05, 
j3 — 1 — Of, and sizes L as specified [ see Eqs. Q and ^ ] . The 
points are from numerical simulations, A'^sam ~ samples. 
The dashed (red) line corresponds to the minimum allowed 
activity (see text). 



Table I. For systems with L sites, and a, j3 as specified, 
(a_t) = {A)/(L + 1) is average current activity [normalized 
as in Eq. (g])]; A = [(ai) - (az,)^]^/^; 'Th' stands for results 
from theory [Eqs. (|29p - (|34[) [. and 'N' for results of numerical 
simulations with A'^sam = 10^, iVset = 10 (see text). 



L 


Type 


(at) 


A 


Skew 






a 


= /3 = 1/2 




64 


Th 


1/4 


0.0305279 . . . 







N 


0.24993(2) 


0.030521(22) 


-0.0008(22) 


128 


Th 


1/4 


0.0218399 . . . 







N 


0.24998(2) 


0.021850(13) 


-0.0005(20) 


256 


Th 


1/4 


0.0155338 . . . 







N 


0.24998(1) 


0.015536(14) 


0.0017(26) 






a = 


1/4, /3 = 3/4 




64 


Th 


3/16 


0.0351324 . . . 


0.01648 . . . 




N 


0.18744(3) 


0.035138(16) 


0.0162(22) 


128 


Th 


3/16 


0.0250771 . . . 


0.01134. . . 




N 


0.18747(3) 


0.025080(17) 


0.0118(22) 


256 


Th 


3/16 


0.0178161... 


0.00790 . . . 




N 


0.18747(2) 


0.017823(13) 


0.0071(23) 


a = /? = 1 


64 


Th 


0.255814.. 


0.0311221 . . . 


0.000091968 . . . 




N 


0.25583(4) 


0.031183(20) 


-0.0017(19) 


128 


Th 


0.252918 . . 


0.0220529 . . . 


0.000016219 . . . 




N 


0.25293(2) 


0.022083(15) 


0.0005(27) 


256 


Th 


0.251462 . . 


0.0156096 . . . 


0.000002864 . . . 




N 


0.25146(2) 


0.015618(8) 


0.0008(12) 




-0.05 0.05 
X 

Figure 4. (Color online) PDF for the reduced variable x = 
{A - {A))/(L + 1) for a = /3 = 1, L = 256 [see Eqs. ©, (p9)l . 
and (|34p [. The points are from numerical simulations, A'^sam = 
10*' samples. The full line is a Gaussian fit to simulation data. 



Table II. For systems with L = 16 sites, a = /3 = 1, re- 
sults for total activity: (a) including {{A")) and (b) not 
including ((A'")) injection and ejection bond contributions. 
A = [{A^) - {Af]^^^; 'Th' stands for results from theory 
[Eqs. and (Ol, and 'N' for resuhs of numerical 

simulations with A'^sam ~ 10^, A'sct ~ 10 (see text). 



Type 




A 


Skew 


(a) {A") 


Th 


4.636363 . . . 


1.042933 . . . 


0.003006 . . . 


N 


4.63655(37) 


1.04286(25) 


0.00295(69) 


(b) {A'") 


Th 


4.090909 . . . 


0.982518... 


0.009706 . . . 


N 


4.09113(34) 


0.98246(12) 


0.00959(34) 



be needed in order to verify the rather small predicted 
values. 

In order to clarify the latter point, we made long runs 
with A^sam — 10^ for L = 16 with a — f3 — 1, evaluating 
both (^") and (A'") of Eqs. (HHl-dni), with the results 
displayed in Table [TTl One can see that the numerical 
results clearly confirm the theoretical prediction for the 
skewness, within error bars. 
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Table III. For systems with PBC, L sites, and p as specified, 
(ol) = {A)/L is average current activity [normalized as in 
Eq. ©]; A = [(ai> - (ai)^]^/^; 'Th' stands for results from 
theory [Eqs. I|5ip - (|53p ]. and 'N' for results of numerical sim- 
ulations with A'sam = 10®, A'set = 10 (see text). 



Table IV. For systems with PBC, L — 16 sites, and p as spec- 
ified, (ul) = {A)/L is average current activity [normalized 
as in Eq. ©]; A = [(ai> - (az,)^]^/^; 'Th' stands for results 
from theory [Eqs. (|5ip - (|53|) [. and 'N' for results of numerical 
simulations with A^'sam = 10^, iVset = 10 (see text). 



L 


Type 




A 


Skew 


p= 1/2 


64 


Th 


0.253968 . . 


0.0312461 . . . 


0.000130123 . . . 




N 


0.25398(2) 


0.031261(20) 


-0.0005(19) 


128 


Th 


0.251969 . . 


0.0220964 . . . 


0.000022272 . . . 




N 


0.25197(2) 


0.022095(11) 


0.0004(17) 


256 


Th 


0.250980 . . 


0.0156249 . . . 


0.000003875 . . . 




N 


0.25099(1) 


0.015629(10) 


0.0001(26) 


p= 1/4 


64 


Th 


0.190476.. 


0.0231771 . . . 


-0.17946... 




N 


0.19048(2) 


0.023178(11) 


-0.1799(13) 


128 


Th 


0.188976.. 


0.0164837... 


-0.12226... 




N 


0.18898(1) 


0.016492(10) 


-0.1237(30) 


256 


Th 


0.188235.. 


0.0116877... 


-0.08487... 




N 


0.18823(1) 


0.011690(7) 


-0.0833(29) 



C. Periodic boundary conditions 

For systems with PBC, the p ^ 1 — p duality exhibited 
in Eqs. ([S^ - fSS]) ensures that all relevant aspects can be 
investigated by considering only, e. g-, p < 1/2. We 
evaluated the first three moments of the activity PDF 
for densities p = 1/2 and 1/4, and varying system sizes 
L. Results are given in Table Hill together with the the- 
oretical predictions summarized in Eqs. (j51|) - (j53p . Sim- 
ilarly to the open boundary case, the agreement is very 
good, except that accuracy for the skewness is compar- 
atively low. For p — 1/2, in particular, comparison of 
simulational data with predicted values of skewness [ see 
Eq. ((55)) I is difficult because the latter are rather small. 

We then made long runs with iVsam — 10^ for L ^ 16 
for both p = 1/2 and 1/4, with the results displayed 
in Table IIVI In this case the numerical results clearly 
confirm the theoretical prediction for the skewness: very 
accurately for p — 1/4, and within reasonable error bars 
for p = 1/2. 

As a consequence of Eq. ([57)) . for p 7^ 1/2 the skewness 
of the PDFs for PBC tends to grow, in absolute value, 
as \p — 1/2| increases. For system sizes L large enough, 
so that the set of (discrete) allowed activity values can 
be reasonably well represented as a continuum, one can 
approximate the respective PDF by a Gaussian with a 
skew-inducing perturbation: 



Pix) 



G{x',<j) , 



(59) 



where 



xq, G{x' , <t) is a Gaussian curve centered 



aX X — xq, with width cr; a, xo, and a are adjustable 
parameters, with the proviso \a\ <C 1. 



Type 

Th 
N 



0.266667.. 
0.26668(1) 



A 

P=l/2 
0.0623610 . . 
0.062350(8) 



Skew 

0.005140 . . . 
0.00515(72) 



1/4 



Th 
N 



1/5 
0.20002(1) 



0.0443203 . . 
0.044313(4) 



-0.455600 . . . 
-0.45578(30) 
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Figure 5. (Color online) PDF for the reduced variable x = 
{A - {A))/L for system with PBC, p = 1/4, L = 256 [see 
Eq. Q [. The points are from numerical simulations, A'sam = 
10" samples. The full (blue) line is the perturbed Gaussian 
distribution given by Eq. (|59p . adjusted to the numerical data. 
The dashed (red) line is a pure Gaussian fit to the same data, 
and is shown for reference. 



In Fig. [S] we show the results of fitting data for L = 
256, p = 1/4 to Eq. (inni)- The adjusted estimate of 
a = -0.0120(25) corresponds to 5 = -0.072(15), which 
compares favorably with the values given in Table lllli 
albeit with a somewhat large uncertainty. The adjusted 
width a = 0.01171(3) is in very good agreement with the 
A values (both from theory and from numerics) given in 
the same Table. 

For a graphic illustration of large skewness with PBC 
(similar to the effect shown in Fig. [3] for open BC with 
a ^ 1, /? sa 1), one needs \p — 1/2| large, as well as 
large L. In the low-density regime, the maximum allowed 
activity is Amax — M , M being the number of particles 
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Figure 6. (Color online) PDFs for the reduced variable 
X = {A — {A))/L for systems with PBC, p = 1/16, and sizes 
L as specified [ see Eq. ((3)l ] . The points are from numerical 
simulations, A'sam = 10® samples. The dashed (red) line cor- 
responds to the maximum allowed activity (apart from small 
corrections of order , see text). 



Figure 7. (Color online) PDF for the reduced variable x = 
{A - {A))/L for system with PBC, p = 1/16, L = 2048 [see 
Eq. (|3]) ]. The points are from numerical simulations, A'sam = 
lO" samples. The full (blue) line is the perturbed Gaussian 
distribution given by Eq. (I59p . adjusted to the numerical data. 



in the system. One has, from Eq. (|?T|) . (A) = M(l — 
p) (1 + C'(i-i)), thus the maximum allowed value of a; = 
{A - {A))/L is p2 (l + C'(i-i)). This is illustrated in 
Fig. [HI for p = 1/16 and various values of (large) L, 
where the sharp cutoff in the forward end of the PDF 
is evident. For the curves shown, the calculated values 
of the skewness are S = -0.58963 . . . , -0.41258 . . . , and 
-0.29023 . . . , respectively for L = 512, 1024, and 2048 
[see Eqs. (ICT])-(IS21) [• The low-density regime with PBC 
is thus the counterpart of the small a, large /3 example 
for open BC of section UlIBI Note that the skewness here 
has the opposite sign to that case. 

We considered the system with p — 1/16, L — 2048 in 
further detail. Fig.[7]shows the results of numerical simu- 
lations, as well as their best fit to Eq. ([5^ . The adjusted 
value of a = —0.043(9) corresponds to S* = —0.27(6), in 
broad agreement with the theoretical prediction. How- 
ever, it can be seen that the perturbative scheme of 
Eq. ([5^ cannot properly account for the long tail at neg- 
ative X. Although we have checked that the tail still 
varies as ~ exp[— 6(x — xq)^], it is not possible to find a 
single set of parameters in Eq. ([5^ which optimizes the 
fit to the latter region, without seriously compromising 
the description of data away from the tail. 

In Ref. [12, an approach similar to the one described 
above was used to calculate current activity PDFs for 
one-dimensional flow of extended objects with exclusion, 
on a ring. Those were compared to results of numer- 
ical simulations of the same quantity, with very good 



agreement (see their Figure 1). Note that their system 
of extended objects corresponds to an effective lattice 
of L = 35 sites with M — 15 point particles. For this 
case Eqs. (l5l])-(IS31) give, e.g., S = -0.015 .... So this is 
a regime where deviations from Gaussianity are rather 
small for activity distributions, which is visually con- 
firmed by the aspect of their figure referred to above. 

IV. DISCUSSION AND CONCLUSIONS 

The investigation of current fluctuations is usually car- 
ried out by examining the total charge crossing a given 
bond, during a long time interval in the steady state 
regime [7Hl3|. Conversely, the snapshot nature of cur- 
rent activity means that samples are collected at a fixed 
time; in this case, the non-trivial features arise when one 
considers the corresponding global quantity, i.e., the sum 
of contributions from all bonds in the system, at a given 
instant. When focusing on fluctuations of either quantity, 
the task for currents is made more involved by the need to 
subtract a background term which grows linearly in time, 
and (depending on the speciflcs of the case), additional 
sublinear terms as well. Current activity fluctuations, on 
the other hand, can be adequately sampled by resorting 
to large enough systems, so that discrete-lattice effects 
are minimized. While, as emphasized earlier, these two 
quantities are distinct in character, both can yield phys- 
ical insights into the properties of flow with exclusion. 

The results of the numerical simulations, reported in 



12 



Section IIIIl provide unequivocal support to the predic- 
tions summarized in Eqs. (0IJ1) - (PT|) [for the cases where 
Eqs. (1221) and §^ apply], and in Eqs. We 
thus collected further information, regarding points of 
the a — (3 phase diagram for open BC not investigated 
in Section lllli by applying the general expression for the 
quantities Zl, Eq. ([22), directly to Eqs. (HHIl-lEIl)- The 
results are as follows. 

We focused on the large-L dependence of the cumu- 
lants Cn 

titles A = C^'^iL 



3, as reflected in the associated quan- 
1) and 5* = Cs/C^^^, as well as on 
the sign of the latter. 

From numerical evaluations of Eqs. (P9)) - (PT|l using 
Eq. ([211) we found A oc L^^/^ everywhere, and strong 
evidence for 



S cc L 



-3/2 



(60) 



in the limit of large L, where km = — 1 for some extended 
regions of the a — /? phase diagram, km — for other re- 
gions, as well as special lines, and km — 1 for some special 
points (see FigurelH). Eq. ((60)) corresponds to C3 oc L^''"^ 
at large L. All this is consistent with results obtained by 
using in Eqs. (I2ni)-(I2I]), instead of Eq. ^ for the 
asyrnptotic approximations to it provided by Refs. 0, [ij 
andlil (see, e.g., Eqs. (8) and (9) of Ref. [H). This 
gives, for the low-current phases [a < 1/2, /? (low den- 
sity) or (3 < 1/2, a (high density)] and the coexistence 
line {a — /3 < 1/2) the results 

C2 = L^(l-3/i) + 0(l) ; 

C3 = L^l{l-4fl){^-5^i) + o{l) , (6i) 

where /i = mm{a, l3) [1 — min(a,/3)]. Thus, in these 
cases, away from the dual lines n — 1/5 the large-L 
skewness is 5 oc L~^^^ and positive (negative) for fi 
less (greater) than 1/5. At = 1/5, corresponding to 

min(a,/3) = ^ ^ "f^ = 0.27639..., the vanishing of 
the 0{L) term in C3 makes l^] < 0{L-^/^). These 
agree with the signs and L-dependences of S from the 
numerical calculations using Eq. (j32l) . gathered in Fig- 
ure [S] There the long-dashed lines are the dual lines 
where fi = 1/5. They cross the a + P = 1 line at 
af3 = 1/5. At this point, the term in C3 contributing 
toS'ocL^^/^ indeed vanishes, as anticipated in Eq. (ITSl) . 
In fact, Eqs. (|5T|) reproduce Eqs. dT5|) and (ITSl) . every- 
where along the a + /S — 1 line, including a — /3 — 1/2 
where, from Section III Al we know that S vanishes. 

For a, P inside the high current phase boundaries, us- 
ing the asymptotic Zl (Eq. (8) of Ref. IB') in Eqs. (|29|) - 
(I^Tj) gives S oc L^'^l'^ ^ in agreement with the results ob- 
tained using the full form, Eq. (I32p . Furthermore, the 
full form shows that near (q;,/3) = (|,^), S > for 
small L, crossing over to negative values as L increases. 
For example, at a — (3 — 0.55, S goes through zero for 
L « 300. The L~^/'^ dependence involves an exception to 
the normal subdominance of injection and ejection con- 
tributions: both the approximate and exact procedures 
show that without the injection and ejection contribution 
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Figure 8. (Color online) The sign of skewness S in the various 
regions of the a — j3 phase diagram is shown. The large-L 
dependence, 15*1 oc L~^, is: in the low-current phases [a < 
1/2, P and j3 < 1/2, a] x = 1/2, except for the long-dashed 
(blue) lines, including their extremes at /? = 1 and a — 1 
[x = 3/2 ] (see text for explanation of these lines) . In the high- 
current phase a, P > 1/2, x = 3/2. Full (red) lines separating 
high- and low-current phases: x = 3/2. Full squares: x = 5/2. 
On the a = 1 and 13 = 1 lines, S has the same sign and L- 
dependence as in the respective adjacent regions, except for 
the points marked by full squares. The circle marks (a, /?) = 
(1/2, 1/2) where 5" = [see Eq. dTTJ]. The short-dashed line 
is the coexistence line between high-and low-density phases. 
Because of particle-hole duality, S is the same for pairs of 
points symmetric with respect to the a = j3 line. 



[i.e., considering the (A'") of Eqs. (H^-lEIl) ], S would 
have been ©(L^^/^); an associated relationship to cumu- 
lants of A' helps to explain the result 5* oc L~^/^ at the 
special boundary point a = (3 = 1 \ see Eq. ([551) ] . 

It is remarkable that the asymptotic approxima- 
tions for Zl are sufBciently accurate, when inserted in 
Eqs. ([^^ - ([HT]) . to give the correct asymptotics for the 
cumulants and skewness, given the high orders of cancel- 
lations in 1/L expansions often occurring. Direct calcu- 
lation of activity cumulants using asymptotic PDFs given 
at the end of Section lllB 3l are less accurate: for example, 
the Laplace approximation mentioned there gives, inside 
the low current phases, S oc L~^/^, but does not give the 
crossover near the long-dashed (blue) lines of Fig. [51 

Correspondence with select results for standard cur- 
rent studies is as follows. In Ref. [13, current statistics 
were considered for open BC. At large L, for the PDF 
of total charge fluctuations, Q — Q{T) — {Q{T)), ac- 
cumulated up to time T they found critical scaling at 
a = 13 = 1/2 with T^/^^'PTiQ) = f{Q/T^^^'), z = 3/2; 
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elsewhere on the a + /3 = 1 Hne, scaling of the PDF was 
found to be Gaussian. Although their results are strictly 
not comparable to ours, their Figure 4 suggests that, e. 
g., for L — 256 critical scaling would only persist up to 
times r < 3, 000. Whether or not some such signature 
would show up in the (fixed-time) statistics of current 
activity, for times other than those used in the present 
work, is an open question. 

Secondly, in Ref. [l3 it was found that, at the special 
point a = /3 = 1, C3 approaches a small, finite (negative) 
value, as L — 00. On the other hand, their C2 oc L~^/^. 
Thus, the standard skewness |20] would become very 
large for large L. Here, in contrast, the skewness of the 
current activity approaches zero as L cxd everywhere. 
Furthermore, it is positive at the above-mentioned spe- 
cial point. This implies that any relationship of current 
fluctuations to those of the activity would involve subtle 
features. 
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Appendix A: Transfer matrix results for a~\- jS = 1 

The evaluation of Eq. ^ using the representation in 
which T is diagonal yields 



(Al) 



where 

M± = M±(A) 
and 



1 ± Vl + 4a/37 (7 EE - 1) (A2) 



a±(A) = [{ae^P+p)^x±+a|3-f] 



+ ;9(e^" - 1)] 



(A3) 



[(2a*± - l)^i±] 

The first three derivatives of ln(e'*'"^) at A = then pro- 
vide the following leading cumulants (exploiting (0) — 
l,/i-(0) = 0): 

Ci = (L + l)a/3 

C2 = La/3 - (3L + 1) {aPf 

C3^Lal3- (9L - 1) {a^f + 4(5L - 1) {a^f . (A4) 

These results generalize Eqs. (fT4)) - (|T6|) and show that 
Eqs. dm), HSl), and ^ apply for L > 1, a/3/(l - 3a/3), 
a/3/(l — 5aP) respectively. 

The exact cumulants given above can alternatively be 
obtained from the forms provided for general a, (3 in 
Section HlB 11 taking Eq. ([55]) into account. 



Appendix B: {A ), (A^) for general a, P 



We first consider the product of bond terms 



L-l 



L-1 



L-l 



A'' = E (1 - E (1 - C™) Cm+l E (1 - Cn) Cn+1 

m— 1 n— 1 

= E ^^"i" ■ (-^1) 



e=i 



(£, m, n) identifies a particular term in the triple sum, 
and corresponds to a point in an (L— 1) x (L — 1) x (L— 1) 
cube in 3-dimensional £mn space. 

We need to distinguish and enumerate the following 
cases (a )-(/ ), of which (a )-(e ) involve repeated site la- 
bels: 

(a) {£mn) = 1 < £ < L - 1: Na = L - 1 points 

(on the diagonal of the cube). Tm = (1 — 0)"^ C|+i — 
(1-0)0+1- 

(6) {£mn) = {£,£,£+!), 1 < £ < L - 2: L-2 
points (on subdiagonal line). Tm+i = (1 — 0)^ Cf+i (1 ^ 
0+1) Ce+2 = . 

This is one example out of 6 lines of the type [£mn) — 
{£,£,£ zL 1) and cyclic variations, containing altogether 
Ni, — 6{N — 2) points each giving zero contribution. 

(c) {£mn) — (£, £±1, £t^), (^ varying; £, to, n all differ- 
ent): with cyclic variations, there are 6 lines containing 
points of this type, altogether making Nc = 6{L — 3) 
points, all with Tf mn — 0. 

(d ) {£m to) : as to vary, this corresponds to a plane 
slicing the cube. It contains (L — 1)^ points of which 
L — 1 are on the line {£mn) — {£,£,£) specified in (a), 
and 2(L - 2) are on lines {£,£ ± 1,£±1) of the type 
specified in (6). So there are (L- 1)^- (L- l + 2(L-2)) 
"new" points associated with the plane. Allowing for 
cyclic variations there are three such planes, containing 
altogether Nd = 3(L — 2){L — 3) "new" points, each with 

Tgrnm = (1 - Ce)Q+l (1 - Cm) Cm+1 {m:^ £, £± 1). 

(e ) {£, TO, TO ± 1) with £, to varying (and cyclic varia- 
tions). There are six planes of this type. Discounting 
points already accounted for, the new planes contain a 
total of Ne = 6(L — 3)(L — 4) "new" points, at each of 
which Timn vanishes. 

(/ ) The remaining points {£ m n) have no pair of coor- 
dinates equal or differing by ±1, so the associated Timn 
has no shared site labels. The total number of such 
points is Nf = (L- 1)^ - {Na + Nb + Nc + Nd + Ne) = 
(L-3)(L-4)(L-5). 

The average of the addition of the non-zero contribu- 
tions to A'^ from (a), (d), (/ ) gives 

{A'') = Na{il-Q)Ci+i) + 

+ Nd ((1 - 0) 0+1 (1 - Cm) Cm+l) + 
+Nf ((1 - Q) Q+l (1 - Cm) Cm+1 (1 - Cn) Cn+l) ■ (B2) 

Here no pair from £mn are equal or differ by ± 1 , so 
we can reduce Eq. (jB2l) using Eqs. (|23l) . and their 
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generalization for this case 



Appendix C: (A^) for PBC 



Zl- 



((l-0)0+l(l-U)U+l(l-Cn)Cn+l) = ^ . (B3) 

We find 

{A'')^{L-3){L-4){L-5)^ + 

+i{L - 2)(L - 3) ^ + (L - 1) ^ . (B4) 

The generahzation for inclusion of injection and ejection 
contributions in the third moment of the activity involves 
similar steps to those for the second moment, given in 
Eq. (HID. We first write 



=Tq + 3Ti + 3r2 + Ta 



where 



T„ = (A')" (aCi+/3(l-CL)) 



3— n 



(B5) 



(B6) 



Using («Ci)" = Ci, («Ci)" (1 - Ci) C2 = (m > 1) 
and similarly for /3 (1 — Cl), one finds 

(To)^(a^ + /3^)%i+3(« + /3)^, 

(Ti) - (L - 2)(a + P) + 2(i - 3) ^^"^ 



(r2)=2(L-2)%i + 2(L-3)(L-4)%^ . (B7) 



Along similar lines to those followed in Appendix |BJ 
we consider the product: 



£—1 tn—l n—1 

(Cl) 

with (i+i = Cl- 

The cases involving site labels are the same set {a)-{f) 
described in Appendix [Bj only the enumerations vary. 
Recalling Eq. (jB2p . one has now, adding up the nonzero 
contributions: 

{A^)=Nl{{l-C,)Q+i) + 

+ 7V,^ ((1-0)0+1(1- Cm) + 
+iVP ((1 - 0) 0+1 (1 - Cm) Cm+1 (1 - Cn) Cn+l) • (C2) 

The numerical coefficients are as follows: 

(a) {(.mn) = (£,^,£), 1 < i < L: NP = L points. 
(d) (imm): = 3L(L - 3) points. 
(/ ) The points (i m n) which have no pair of coordinates 
equal or differing by ±1: — L{L — 4)(i — 5) points. 

One has: 



((1-0)0+1) = 



lCm 



(C3) 



((1 - 0) 0+1(1 - Cm) Cm + l) = ^^1^ ; (C4) 



With (Ts) = {A") from Eq. (|B4)) . 
(A3)^(L-l)(L-2)(L-3)^ + 

+3(L - 1)(L - 2 + a + /?) + 

2 , o2\ ^L-1 



+ (L-l + a'+r) 



(B8) 



((1 - 0) Cf+1 (1 - Cm) Cm+1 (1 - Cn) Cn+l) — 

So, recalling Eq. 

(A^) ==ao(l + 3ai + aia2) . 



(C5) 
(C6) 
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